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Abstract 

Using the theory of negative association for measures and the notion of random weak limits of 
sparse graphs, we establish the validity of the cavity method for counting spanning subgraphs subject 
to local constraints in asymptotically tree-like graphs. Specifically, the corresponding free entropy 
density is shown to converge along any sequence of graphs whose random weak limit is a tree, and 
the limit is directly expressed in terms of the unique solution to a limiting cavity equation. On a 
Galton- Watson tree, the latter simplifies into a recursive distributional equation which can be solved 
explicitely. As an illustration, we provide an explicit-limit formula for the 6— matching number of an 
Erdos-Renyi random graph with fixed average degree and diverging size, for any 6 £ N. 



1 Introduction 

The general framework we consider is that of a finite graph G = {V, E), in which spanning subgraphs are 
weighted according to their local aspect around each vertex as follows : 



M^) = Y[t^.{FnE,). (1) 

iev 

Here, a spanning subgraph (V, F) is identified with its egde-set F <Z E, and each fii is a given non-negative 
function over the subsets of Ei {e G E; e is incident to i}. We call fj, the global measure induced by 
the local measures fii,i G V. Of particular interest in combinatorial optimization is the number 

M(G) = max{|F| : F e supp(Af)}, (2) 

which is the maximum possible size of a spanning subgraph F satisfying the local contraint ^i{FC\Ei) > 
at every node i €z V . More generally, counting the weighted number of spanning subgraphs of each given 
size in G, i.e. determining the generating polynomial 

ZiG;t) = (3) 

FCE 
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is a fundamental task, of which many combinatorial problems are special instances. Intimately related 
to this is the study of a random spanning subgraph J- sampled from the Gibbs-Boltzmann law : 



Z{G;t) ' 



(4) 



where i > is a variable parameter called the activity. In particular, the expected size of J- is called the 
energy U{G;t) and is connected to Z{G;t) via the elementary identity 



Our concern is the behavior of these quantities in the infinite volume limit : \V\ oo, \E\ = 0(|y|). 

Originating from spin glass theory |26| . the cavity method is a powerful nonrigorous technique for 
evaluating such asymptotics on graphs that are locally tree-like. Essentially, the heuristic consists in 
neglecting cycles in order to obtain an approximate local fixed point equation for the marginals of the 
Gibbs-Boltzmann law. Despite its remarkable practical efficiency and the mathematical confirmation of 
its analytical predictions for various important models (JU [1] [501 UHl IZl 113 US] > this ansatz is still far from 
being completely understood, and the exact conditions for its validity remain unknown. More precisely, 
two crucial questions arise in presence of cycles : 

1. convergence : is there a unique, globally attractive fixed point to the cavity equation ? 

2. correctness : if yes, does it have any relation to the Gibbs-Boltzmann marginals ? 

In this paper, we exhibit a general condition under which the cavity method is valid for counting 
spanning subgraphs subject to local constaints. Specifically, we positively answer question 1 for arbitrary 
finite graphs (Theorem [2|) , under the only assumption that each local measure enjoys a certain form of 
negative association which we call the cavity-monotone property, and which simply boils down to ultra-log- 
concavity in the exchangeable case. Regarding question 2, we use the framework of local weak convergence 
[5J 13] and the notion of unimodularity [3] to prove asymptotical correctness for any sequence of graphs 
whose random weak limit is concentrated on trees (Theorem|4]). This includes many classical sequences of 
diluted graphs, such as random d— regular graphs, Erdos-Renyi random graphs with fixed average degree, 
or more generally random graphs with a prescribed degree distribution. In all these examples, the limit 
is a unimodular Galton- Watson (UGW) tree. Thanks to the distributional self-similarity of such a tree, 
the cavity equation simplifies into a recursive distributional equation which may be solved explicitely. As 
a motivation, let us first describe the implications of our work in the special case of b-matchings. 

A famous combinatorial structure that fits in the above framework is obtained by fixing 6 £ N and 
taking ^i(F) = 1(|F| < &) for alH £ F : the induced global measure fi is then nothing but the counting 
measure for 6— matchings in G, i.e. spanning subgraphs with maximum degree at most b. The reader 
is refered to the monograph |30| for a complete survey on &— matchings, and to [24) for the important 
case of matchings (b — 1). The associated quantities Mi,(G) and Zi,(G;t) are important graph invariants 
respectively known as the b— matching number and b— matching polynomial. Determining Zi(G;i) is a 
classical example of a computationally hard problem j32i . although efficient approximation algorithms 
have been designed [SI [S]. The mathematical properties of Zb{G;t) have been investigated in detail, 
notably in the case 6 = 1 for the purpose of understanding monomer-dimer systems [21ill0| . Interestingly, 
the geometry of the complex zeros of Zb{G, t) has been proven to be quite remarkable (see [21] for 6 = 1, 
[29] for b — 2, and [34] for the general case). Regarding Mi{G), the first results in the infinite volume 



U{G;t) = t^\ogZ{G;t). 



(5) 



2 



limit were obtained by Karp and Sipser [23' for the Erdos-Renyi random graph G„ with average degree 
c > on n vertices : 

-Mi(G„)-^l-^ ^ , (6) 

where ic G (0, 1) is the smallest root of i = e^"^"^ ^* . The analysis has then been extended to random 
graphs with a log-concave degree profile [I^ . and finally to any graph sequence that converges in the 
local weak sense [T^. Contrastingly, only little is known for b > 2 : to the best of our knowledge, the 
limit of Y^Mi,{Gn) is only known to exist in the Erdos-Renyi case [TH], and could not be explicitely 
determined. As a special case of our main result, it will follow that 

Theorem 1 (6— matchings in locally tree-like graphs) For any sequence of finite graphs {Gn)neN 
satisfying |i?„| — 0(|V^i|) and whose random weak limit C is concentrated on trees, the limits 

fi,{C) -.^ lim logZfc(G„;i) and m.b{C) -.^ lim 

n-^oo |V„| n-i-oo | V/„ | 

exist and depend only on the random weak limit C When C is a UGW tree, we have the explicit formula 

mb{C) min jfe- ^gtis) ~ ^{gbofb)is) + ^fb{s){fb o /h)(s)| , 
sG[04] 1^ 2 2 2 ) 

where c, /, g are defined in terms of the degree generating function ip{s) = tt^s'^ as follows : 
c = cj,{l), Ms) ^-2^ ^, and gbis)^}^ . 

A;=0 ■ fc=0 

Moreover, any s where the above minimum is achieved must be a root of s — (fb o fb){s). 

For example, in the case of Erdos-Renyi random graphs with average degree c > on n vertices, the 
random weak limit C is a.s. the law of a UGW tree with Poisson(c) degree distribution, and hence. 



1 . , X a.s. . ( , b , , b 



c 



-Mb{Gn) niin <^ b- -Vb+i{s) - 7;{Vb+i ° Vb){s) + -ipb{s){ipb o (pb){s) } , 

n n->oo se[o,i] [_ 2 2 2 ' 



where we have set 

fc-i I, 
cs'' 



k=0 



Since any s where the minimum is achieved must satisfy s = (ipb ° ipb){s), we recover exactly ^ in the 
special case of matchings {b — 1). 

The paper is organized as follows : in section [21 we recall the basic notions and properties pertaining 
to measures over subsets, which will be of constant use throughout the paper. In section[31 we define and 
study the cavity equation associated to a finite network. In section |4l we extend the results to infinite 
networks that arise as local weak limits of finite networks. Finally, section[5]is devoted to the study of the 
cavity equation in the limit of infinite activity, and to its explicit resolution in the case of 6— matchings. 
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2 Preliminaries 



In this section, we define the important notions pertaining to (non-negative) measures fj, over the subsets 
of an arbitrary finite ground set E. Later on, these will be specialized to the local measures {iJ,i)i^v 
attached to the vertices of a graph G. First, /i is caracterized by its multivariate generating polynomial 

Z(w) = ^m(F)w^, (7) 

F<ZE 

where w = {we)eeE and = J^gg^ Wg. Since Z is affine in each w^, e E E, it can be decomposed as 

Z(w) = u.eZ/'=(w') + ^\"(w'), (8) 

where w' = {wf)f^e and Z^*^, Z/*^ are the multi- affine polynomials with ground set E\e respectively 
obtained from Z by setting the variable We to (deletion) and differentiating with respect to We (con- 
traction). By definition, the cavity ratio of the pair (/i, e) is then simply the multi-affine rational function 

r-(w') = V^-T- (9) 

When positive values are assigned to the variables (a so-called external field), we may consider the 
probability distribution PJ[(T = F) = fi{F)w^ / Z (w) . A quantity of interest is the expected size of J- 
when viewed as a function of the external field. We call this the energy : 

U,{xv) = E-[|^|]. (10) 

From the decomposition ([5]), it follows immediately that 

p-(ee.F) = ^ / M > hence C/.(w) - V ^ / m ■ (H) 

l + w;er^(w')' '^^ ^^l + w,T-^{^') 

Note that the supremum of the energy is exactly the rank of fi : lank(fi) = max {[F|; F S supp (fi)} . The 
following properties will be of crucial importance throughout the paper. 

Definition 1 (Cavity-monotone measures) The measure fi is called 

• Rayleigh if every two distinct ground elements e ^ / are negatively correlated in J- : 

Vw > 0, P; (e G ^, / e ^) < P; (e e F) P; (f^F). 

• Size-increasing if every ground element e positively influences the total size \F\ : 

Vw > 0,E; [\F\^eeT)] > E; [|^|]P;(e e F). 

• Cavity-monotone if its satisfies fi(9) > and the two above properties. 

Rayleigh measures were introduced in the context of matroid theory [33j , but soon found their place in 
the modern theory of negative dependence for probability measures [27\ I22j . Cavity- monotone measures 
will play a major role in our study, for the following elementary reason. 
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Lemma 1 (Monotony of energy and cavity ratios) 

/u(0) > •<==^ the cavity ratios F^, e d E are well-defined on [0, cxd) x . . . x [0, cxd). 

/i is Rayleigh the cavity ratios rjj,e £ i? are non-increasing in each variable. 

fi is size-increasing the energy is increasing in each variable] 

for each e £ E and w > 0,t ^ trj^(iw') is increasing. 

Proof. Differentiating tlie corresponding quantities and playing with the definition of easily yields 



dWe 



dtV^ (tw') 

'^Q^ = K [mT-(eer)] - m] p;r(e e ^)) /^^eP;: ^ -^)' . 

Remark 1 (Matroids) Interestingly, the .support of a cavity-monotone measure admits a remarkable 
structure : it follows from \39\. Theorem 4-.6] that for fi Rayleigh with fi{9) > 0, I — supp fi is a matroid: 



• I is not empty ; 

• If B e I and AC B, then Ae I ; 

• If A,B e I and \ A\ < \B\, then 3e <E B\A such that AUe e I. 

The cavity-monotone property admits a particularly simple caracterization in the important case 
where /i is exchangeable, i.e. fJ,{F) — c(\F\) for some non-negative coefficients c(0), . . . ,c{m),m ~ \E\: 

Lemma 2 (The exchangeable case) An exchangeable measure fi is cavity-monotone if and only if 

1. c is log-concave, i.e. c^{k) > c{k — l)c(fc + 1) for all < k < m, and 

2. the support {0 < k < m : c{k) > 0} is an interval containing and 1. 

In particular, so is the measure fJ-{F) — l(|^|<b) describing the local constraints of a 6— matching. 

Proof of Lemma\^ The result essentially follows from the work of Pemantle 27!. Indeed, Theorem 
2.7 therein guarantees that is Rayleigh if and only if the sequence c is log-concave and its support 
is an interval. That the latter must contain is nothing but the last property in the definition of a 
cavity-monotone measure. That it is not reduced to is imposed by the strict inequality in the size- 
increasing property. Conversely, let us show that any exchangeable measure /x with c(0) > and c(l) > 
is indeed size-increasing. Fix an external field w > 0. By Lemma 2.9 in [27], the law obtained from P^ 
by conditionning on the event {|J^| = fc} is stochastically increasing in k. By Proposition 1.2 in [27], this 
implies in particular that for every e C E, the following weak inequality holds : 

E; >E;[|^|]. 

Note that the condition c(l) > guarantees that this conditional expectation is well-defined. Since we 
have not yet used the fact that c(0) > 0, the above inequality remains true if one changes the coefficient 
c(0) to 0. Setting it then back to its initial (positive) value does not affect the left-hand side, but strictly 
decreases the right-hand side, hence the desired strict inequality. □ 
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3 The cavity equation on finite networks 



Let G — (V, E) be a finite graph at tlie vertices of wlrich some local measures iJ.i,i £ V are specified. We 
call the resulting object a network. A configuration x is an assignment of numbers Xi^j > to every 
oriented edge i ^ j £ E. Starting from a configuration x, we define a new configuration y = rG'(x) by 

y,^, ^r%ixk^,:kedi\j), (12) 

where di denotes the set of all neighbors of i. Each Xi^j may be thought of as a message sent by i to j 
along the edge ij, and Tq as a local rule for propagating messages. For t > 0, the fixed point equation 

x-trG(x) (13) 

is called the cavity equation at activity t on the network G. Its relation to the global measure fj, induced 
by the {pi)i^v is revealed by the following well-known result. 

Lemma 3 (Validity on trees) Assume that G is finite and acyclic. Then, for every activity t > 0, 

1. convergence : the cavity equation admits a unique solution x(t), which can be reached from any 
initial configuration by iterating tVc a number of times equal to the diameter of G ; 

2. correctness : for every i G V , the exact marginal law of J- H Ei under the Gibbs-Boltzmann law 
¥q is given by directly imposing the external field {xj^i{t) : j £ di} onto the local measure iii. 

The important consequence is that on trees, the energy U (G; t) can be determined using only local 
operations : 

where the second equality is obtained by applying (1111) to each iii,i £ V. 

Proof of Lemma . When i is a leaf, the message yi^j defined by equation (|12p does not depend at 
all on the initial configuration x. Iterating this argument immediately proves the convergence part, and 
we now focus on correctness. Let G = (V, E) be a finite tree, o a vertex, and i a neighbour of o. We let 
Gi^o denote the subtree induced by o and all vertices that the edge io separates from o. Now assume 
that G is equipped with local measures, and let Gi^o inherit from these local measures, except for /io 
which we replace by the trivial local measure with constant value 1. With these notations, any spanning 
subgraph F C E can be uniquely decomposed as the disjoint union of a subset I C Eo and a spanning 
subgraph Fi on each Gi^o, i £ 9o, with io ^ Fi. Thus, writing for the global measure on the network 
G, we have for any t > 0, 

f,G{F)t^^^ = /^o(/)n*"'''^VG.^„(i^.Uzo)J]tl^'l;,G.^„(F,). 

Fixing / and summing over all possible values for Fi,i £ 9o, we obtain 

P*.(-Fn£;o = /) = G/io(/)nPG._„(»°e-^)n^G.^o(*°^-^)- 

rot 



^Vo(/)np* 
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where C, C" are normalizing constants that do not depend on /. This aheady proves that the law of 
J- n Eo can be obtained from the local measure /io by imposing on each edge io e Eo the external field 

P^. (io e J") 

In turn, this ratio can now be computed by applying the result to the vertex i in the network Gi^o '■ 

which shows that the configuration x(t) defined on G by satisfies the cavity equation (|I3p . □ 

There are two distinct parts in Lemma [3] : convergence and correctness. As we will now show, the 
former extends to arbitrary graphs under the only assumption that each local measure is cavity-monotone. 
Henceforth, such a network will be called a cavity-monotone network. 

Theorem 2 (Convergence on finite cavity- monotone networks) On a finite cavity-monotone net- 
work, the cavity equation admits a unique, globally attractive fixed point x{t) at any activity t > 0. 

Proof. Fixing t > and starting with the minimal configuration x*^ := 0, we set inductively 

x'^+i(t) :=trG(x''^W), 

for all fc e N. By Lemma (TJ the Rayleigh property of the local measures fj,i,i G V ensures that Tq is 
coordinate- wise non-increasing on the space of configurations. Therefore, the limiting configuration 

x-(t) lim tx2'=(i) and x+(t) := lim J, x^'^+i (t) (16) 

exist, and any fixed point x — tTci^) must satisfy x~(t) < x < x+(t). Moreover, Tq is clearly continuous 
with respect to the product topology on configurations, so that tTG{x~{t)) — x+(t) and trG(x+(t)) = 
x^(t). Thus, the existence of unique globally attractive solution to ([T3|) boils down to the equality 

x-(i)-x+(t). (17) 

Now applying (fTTT) to the local measure at a fixed vertex i G V yields 

Summing over all vertices i gV , we therefore obtain 

x^^iit)xt_^j{t) x+^^xT_^^{t) \ 



This implies (jl7p , since by Lemma [1] each U^. , z £ V" is strictly increasing in every coordinate. □ 
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4 The limit of infinite volume 



In the previous section, we have established existence and uniqueness of a cavity solution on any finite 
cavity-monotone network. Our concern now is its asymptotical meaning as the size of the underlying graph 
tends to infinity. Following the principles of the objective method [4], we will replace the asymptotical 
analysis of our finite networks by the direct study of their infinite limits. 



4.1 Random weak limits 

We first briefiy recall the framework of local convergence, introduced by Benjamini and Schramm [8] and 
developped further by Aldous and Steele 4 . Examples of successful uses include jH [9j [HI EH HH [TH 
1131 116) . Here, a network will be simply a denumerable graph G = {V, E) whose vertices are equipped 
with local measures ^i,i dV. A rooted network (G, o) is a network together with the specification of a 
particular vertex o Q V, called the root. For e > 0, we write (G",o') = (G, o) if there exists a bijection 
7 : V —i' V that preserves 

• the root : 7(0) = o' ; 

• the adjacency : ij ^ E <;=4> 7(i)7(j) G E' ; 

• the support of the local measures : fJ.i{F) > t^'y(^i-j{j{Fj) > 0, with 7(i^) = {j{i)'y{j) ■ ij £ F}. 

• the values of the local measures, up to e : \fj.'^^^^{'-f{F)) — Hi{F)\ < e. 

We let C/, denote the set of all locally finite connected rooted networks considered up to the isomorphism 

relation =. In the space C?*, a sequence {(G„,o„);n e N} converges locally to (G, o) if for every radius 
e N and every e > 0, there is nk,e & ^ such that 

n>nk => [G„, o„]fe = [G, 0]^., 

where [G, o]^. denotes the finite rooted network obtained by keeping only the vertices lying at graph- 
distance at most k from o. It is not hard to construct a distance which metrizes this notion of convergence 
and turns G into a complete separable metric space. We can thus import the usual machinery of weak 
convergence of probability measures on Polish spaces (see e.g. [IT). 

Uniform rooting is a natural procedure for turning a finite deterministic network G into a random 
element of 5* : one simply chooses uniformly at random a vertex o to be the root, and restrains G to 
the connected component of o. If (G„)„gN is a sequence of finite networks and if the sequence of their 
laws under uniform rooting admits a weak limit C G Vi^G*), we call C the random weak limit of the 
sequence (G„)niEN- In [3]) it was shown that any such limit enjoys a remarkable invariance property 
known as unimodularity : let G** denote the space of locally finite connected networks with an ordered 
pair of distinguished adjacent vertices (G, o, z), taken up to the natural isomorphism relation and endowed 
with the natural topology. A measure C G V^G*) is called unimodular if it satisfies the Mass-Transport 
Principle : for any Borel function /: G** [0, 00], 



C 



.iedo 



(18) 



where we have written C[-] for the expectation with respect to C This is a deep and powerful notion, 
which we will now use to extend the results of section [3] to the infinite setting. 
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4.2 Main result : validity of the cavity method on unimodular trees 



The definition of remains valid for any locally finite network G. When the latter is cavity-monotone, 
the configurations x~(f),x+(t) introduced in the proof of Theorem [5] remain perfectly well-defined, and 
the convergence of the cavity method again boils down to the identity x~(i) = x+(t). However, the proof 
of the latter involves a summation over all edges, which is no longer valid in the infinite setting. Instead, 
the desired x^ it) — x+ {t) will be derived from unimodularity, and will thus hold for any random weak 
limit of finite networks. Indeed, applying the Mass- Transport Principle to the function 

o, I) .- 



yields C \U ^^{x~_^^{t) : i G 9o)] — C (a:^^ (t) : i G 9o)] (/ is Borel as the pointwise limit of continu- 
ous functions). Under the assumption H [rank (/io)] < oo, this expectation is finite, and the size-increasing 
property of /io then implies that £— almost surely, x~^^{t) = xf_^^{t) for all i G do. This automatically 
extends to every oriented edge since under unimodularity, everything shows up at the root (another fruitful 
application of the Mass- Transport-Principle, see |S1 Lemma 2.3]). We state this as a Theorem. 

Theorem 3 (Convergence of the cavity method on unimodular netv^rorks) Let C be a unimod- 
ular probability measure supported by cavity-monotone networks. // £ [rank (/io)] < oo, then the cavity 
equation admits C—a.-s. a unique, globally attractive solution x(t) at any activity t > 0. 



By analogy with formula (jl4l) in the finite case, the (now well-defined) quantity 

(19) 



u{C;t)^^C 



Xi-io{t)Xo^i{t) 



t -\- Xi^o[t)Xo^i{t) 



appears as a natural candidate for the limiting energy of any sequence of finite networks whose random 
weak limit is C Our second result is the validity of this cavity ansatz when C is concentrated on trees. 

Theorem 4 (Asymptotical correctness of the cavity method) Let (G„)„gN be a sequence of finite 
cavity-monotone networks admitting a random weak limit C which is concentrated on cavity-monotone 
trees. Assume that the local rank at a uniformly chosen vertex is uniformly integrable as n ^ oo. Then, 

> u{C;t). (20) 

\Vn\ n^ca 

If \En \ = 0{\Vn\) and all the local measures take values in {0} U K for a fixed compact K C (0, oo), then 
^logZ(G„;i) C[\og^iom+ f "^^^ds, (21) 

M{Gn) 







^ m{C) := lim t u{C;t). (22) 



\Vn\ 

Remark 2 (Large deviation principle) Integrating \20\) immediately implies that 

J_ ZiGn;t) f' u[C-s) ^ 

|K|^^Z(G„;1) ^ A s 

It will later be checked that t i— > u{C\t) is continuous on M_|_ (Remark\B^. Therefore, denoting by Tn a 
random spanning subgraph with the Gibbs-Boltzmann IuwVq^, Gartner-Ellis Theorem flT^ guarantees that 
obeys a large deviation principle with rate \Vn\ and good rate function y i-^ (y — u{C; e^))~^ ds. 
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4.3 Proof of the main result 



Lemma 4 (Tree approximation) Let (G, o) be a finite rooted cavity-monotone network, and let fc G N. 

// [G, o]2k+2 is a tree, then for every activity t > Q, 

Proof. The proof makes use of a classical ingredient known as the spatial Markov property, which we 
first briefiy recall. Let G = (V, E) be a finite network and let S be an induced subgraph. We let dS denote 
the boundary of S, i.e. the set of edges having one end-point in S and one in S"^. Any boundary condition 
B C dS can be used to assign local measures to the vertices of S, namely fJ.f{F) := /i^ {F U {B D E.i)). 
Note that these local measures differ from the original ones only for vertices that are adjacent to the 
boundary. The resulting network is denoted by S\B. Now, a spanning subgraph F C E is clearly the 
disjoint union of a spanning subgraph Fint of S, a boundary condition B C dS and a spanning subgraph 
i^ext in S'^. The product form of /ic immediately yields : 

P*g(-F = F)= P^|b(-F = Fi„t)P^(-F ndS = B)P*5.|b(-F =. Fe,t). (23) 

In other words, conditionally on the boundary B := H dS, the restrictions of to and S"^ are 
independent with law P^|g and P^cig, respectively. Applying this to the tree S — [G, o]2/c+2, 

E*G[|-Fn£;o|] = ^Gi^<^9S 

BCdS 
BCdS 

where we have applied Lemma [3] to the tree S\B, writing x'-'^^(t) for the unique solution to the cavity 
equation at activity t thereon. But by monotony of the cavity operator, each xl^\^{t), i & do must satisfy 
a;f^o(^) ^ x\^^{t) < xf'!^J'{t). Using the size-increasing property of ^o, we see that 

U,. {^l^oit) ■.^e^o)< U,^ [xi^iit) : . e ao) < U,^ (x^i+/(t) : z G 9o) , 

and re-injecting this into the above equation finally yields the desired inequalities. □ 

Let us now see how Lemma [4] implies the convergence ([20|). Let (G„)„gN be a sequence of finite 
cavity-monotone networks admitting a random weak limit C which is concentrated on cavity-monotone 
trees. Denote by £„ e P{Q*) the law under uniform rooting of G„, so that £„ =^ C. We will use the 
short-hand Uk{G,o) = U^^ {x^^^{t) : i G 9o), and Xk{G,o) for the indicator function that [G,o]2k+2 is a 
tree. Lemma [4] guarantees that for any finite cavity- monotone network G and any vertex o, 

Xfe(G, o)u2k{G, o) < E*c [|^| n ^o] < Xk{G, o)u2k+i{G, o) + (1 - Xk{G, o))rank(/io). (24) 

As functions of (G, o), the left-hand side and right-hand side are continuous on Q^,, since they depend 
only on [G, o]2A;+2- Moreover, both are dominated by (G, o) i— rank(/io) which is assumed to be uni- 
formly integrable with respect to the sequence (/I„)ngN- Thus, their expectation under Cn tends to their 
expectation under >C as n oo. But Xk is zero on the support of C, so we are simply left with 

\C [U,^ {xi%M ■■ ^ e do)] < hminf ^^i^ < hmsup ^^i^ < U [U,^ {xf_t;{t) : ^ e do)] . 

Z n— >oo \Vn\ n-i-oo Vn ^ 



^B)El^s[\TnE.\] 
^B)U,^ (x^^_}^{t):iedo), 
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Since the random weak limit £ is unimodular, Theorem |3] finally implies that both the lower and upper 
bounds tend to C [t/^^, {xi^o{t) : i £ do)] = u{C; as fc — > oo. Note that the requirement i2[rank(/io)] < oo 
in Theorem [3] is here automatically fuUfilled, by the uniform integrability assumption. 

It now remains to show (l?T|) and (P^ . The identity ([5]) implies that for any activity t > and any 
finite network G satisfying /i(0) > 0, 

M '^""'^ ^ M £ ^ r 

Now take G = G„ and let rt oo : the compactness assumption guarantees that log/io(0) is bounded 
uniformly in n, so the first term converges to £[log/io(0)]. As per the second one, it tends to -^^^^^ds 
because of (|20p. provided the uniform domination holds in Lebesgue's dominated convergence Theorem. 
The latter fact is ensured by the first inequality in Lemma [5] below, combined with the compactness 
assumption and the fact that |i?„| = 0(|Kj|). The second inequality in Lemma [5] easily guarantees ((22)) . 



Lemma 5 (Uniform controls for the energy) LetG be a finite cavity-monotone network. As a func- 
tion of the activity t, the energy U{G; t) increases from to M (G). Furthermore, the rate of convergence 
to these two extrema can be precisely controlled : 

vt>o, u{G;t) < tJ2Mp^)MpJ); (25) 

ijeE 

yt>l, U{G;t) > Af(G)--^ (|i?|log2 + Vlog^(M,) ) ■ (26) 
where A(/i) — , with max fi = max{^{F) : F S supp (/i)} and min /i = min{fj,{F) : F £ supp (/i)}. 



Proof of Lemma That the energy increases with the activity is equivalent by ([5]) to the convexity 
oi 9 1-^ log2'(G; e^), a direct consequence of Holder's inequality. This also implies that for any t > 1 

£/(ft<)iog* > iog||i|. 

Clearly Z{G;t) > t^'^^^l min fi and Z(G; 1) < 2l^lmax^ and A(/i) < UievMlJ-i), so & follows. Re- 
garding (1^ . we have for any t > 

U{G;t) = ^P*(*je^) 



ijeE 



< 



< 



^l{F u ij) 



t y max 



where the third line uses the standard inequality < max(|, ^) for any a,b,c,d > 0, and the crucial 
fact that fi{F U ij) > =^ fi{F) > (see Remark [y. 
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5 The limit of infinite activity 



The goal of this final section is to describe the important quantity m[C) introduced in Theorem 2] directly 
in terms of a certain local equation which we naturally call the cavity equation at infinite activity. This 
will then be used to establish the explicit formulae that have been mentioned in the introduction. 

5.1 The cavity equation at infinite activity 

Let G be a cavity-monotone network. From Lemma [1] it follows that (t,x) H- tVQ{tii) is increasing in t 
and decreasing in x. We may thus define a limiting cavity-operator by 

rG(x) := lim ttrG(tx). 

t—^oo 

By monotony, Vq ■ [0, oo)^ (0, oo]-^ is well-defined without any ambiguity regarding the order in which 
the limits t oo and x — !■ are taken. Note also that Tq can be composed with Tg '■ (0, oo]^ — ?> [0, oo)^, 
yielding a two-step local update rule on (0, oo]^ which will now play a crucial role. 

Theorem 5 (The cavity equation at infinite activity) Let G he a cavity-monotone network on which 
the cavity equation at activity t admits a unique globally attractive fixed point x(i), for every t > 0. Then, 

X := lim t x(t) (27) 

exists in (0,oo]^j and is the smallest solution to the so-called cavity equation at infinite activity on G : 

X- (rGorG)(x). (28) 

In particular, for any unimodular probability measure C concentrated on cavity-monotone networks, and 
satisyfing C [rank(/Xo)] < oo, we have 

1 



m{C) = -C [Uf,^ [Xi^o : i € do)] 



Proof of Theorem\^ By assumption x'^(/;) — s> x(t) for any t > 0, where x'' = and for all A; G N, 

T^+\t)^tTG{t^Y'^ or equivalently, ^^^^^ = Tg (x'=(i)) . (29) 

But (i, x) I— >■ tVQ{tyL) is increasing in t and decreasing in x, so an immediate induction over k shows that 
t I— ;> t~^x'"'(<) and t H> x'"'(i) are respectively non-increasing and non-decreasing, for every fc g N. Thus, 
1 1— > x(i) is non-decreasing, hence the existence of (l27l) . The identity (|28p is then obtained by passing to 
the hmits in ([29]). Finally, if y G (0, oo]^ satisfies y = (Fg o FG)(y), then for every fc £ N and i > 0, 

x''=(t)<y. (30) 

Indeed, the property is trivial when fc = 0, and is preserved from fc to fc -I- 1 because iFG(iw) < Fg(w) 
for any w £ [0, oo)^. Letting k ^ oo and then t ^ oo in ([50]) yields x < y, as desired. □ 

Remark 3 (Continuity with respect to the activity) Incidentally, we have just obtained that t H> 
t^^x{t) and t I— > x(t) are respectively non-increasing and non- decreasing, so that 

< s <t =^ jx{t) < x(s) < x{t). 

This guarantees the continuity of t i-^ x;(t), and hence that of u{C;t), as promised in Remark\^ 
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5.2 Explicit resolution for b— matchings on GW trees 



Many classical sequences of diluted random graphs admit almost surely a particularly simple random weak 
limit £, namely a unimodular Galton- Watson (UGW) tree (see Example 1.1 in [3]). This random rooted 
tree is parametrized by a probability distribution tt G 'P(N) with finite mean, called its degree distribution. 
It is obtained by a Galton- Watson branching process where the root has offspring distribution tt and all 
other genitors have the size-biased offspring distribution tt £ defined by 

\fn e N,7r„ = (n + l)7r„+i/^ fcTTfe. 

k 

Thanks to the markovian nature of the branching process, the cavity equation at infinite activity simplifies 
into a recursive distributional equation (RDE) (see [5] for a survey). Let us describe it and solve it in the 
case of &— matchings, where 6 > 1 is fixed. The local cavity and energy ratios are simply 

. X]/C[Arl:|/|<fc-l riiG/ n rr/ \ S/C[Arl:|/|<f) Hie/ ^« 

r(a;i, . . . = — and U{xi, . . . , Xn) - 



J2lClN]:\I\<b Tliel '^I<Z[N]:\I\<b Ilie/ 

In the infinite activity limit, the local cavity ratio becomes 

Y=;/ N _ S/C[A']:|/|=6-1 Die/ 

L\Xi,...,Xn) — — =; ^ : 

l^I<Z[N\:\I\=b\.\.i£l 

where all conventions regarding degenerate cases are obtained by taking the corresponding limits. Given 
Q e P((0,oo]), we let Q{Q) e 7'([0, 1])) denote the law of T (Yi, . . . , Y^) , where N has law tt and 
Yi, ^2, ■ • ■ are i.i.d. with law Q, independent of N. Similarly, given P € 7'([0, 1]), we let e(P) £ P ((0, oo]) 
denote the law of F (Xi, . . . ,X^^, where N has law tt and Xi,X2, ■ ■ ■ are i.i.d. with law P, independent of 

TV. Thanks to the markovian nature of the GW branching process, the law Q G 'P((0, oo]) of a message sent 
towards the root in the configuration x must satisfy the RDE Q = {Q o Q){Q). Equivalently, P = Q{Q) 
must satisfy P — {Q o 0)(P). More precisely, letting M{P) denote the expectation of (Fi, . . . , Iat), 
where N has law tt and Yi,Y2, . . . are i.i.d. with law Q = 6(P), independent of N, we may rephrase 
Theorem [5] as follows (distributions are here endowed with the usual stochastic order) : 

Corollary 1 to(/:) M{P), where P G P([0, 1]) is the smallest solution to the RDE P = (9 o e)(P). 

The fixed points of o O turn out to be in one-to-one correspondance with the historical minima of a 
certain function H : [0,1] defined in terms of the degree generating function (/)(s) = J^k ^^s^ '■ 

n(s) = k ' !<;(«) -^(s" /)(») + " /)(«), 



c^ k\ ' ^ k\ 

k=0 k=0 



A historical minima of H is a number s <E [0, 1] satisfying H'{s) = and H{t) > H{s) for all t G [0, 



s 



Theorem 6 (Resolution of the RDE) Let si < . . . < Sr denote the historical minima of the function 
H . Then, the distributional equation P = (0oO)(P) admits exactly r solutions, and they are stochastically 
ordered : Pi < . . . < Pr- Moreover, for each 1 <i <r, we have M{Pi) — H{si). 

In particular, m[C) — miniJ, which is exactly the formula given in Theorem[T] Theorem[6]was established 
in [M] for the special case 6=1, but the proof can easily be adapted to the general case. For the sake of 
completeness, we have included a general proof in the Appendix. 
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Appendix : resolution of the RDE 



First observe that the raappmgs 9,8 and M are aU decreasing with respect to stochastic order, and 
continuous with respect to the topology of weak convergence. Note also that cf'{t)t — bf{t), so that 
Wit) = cf'it) ((/ o /)(<) - t) . Thus, Wit) = if and only if (/ o /)(<) = t. 

Lemma 6 (Properties of F, F and U) 

1. Let ixi,...,Xn) e [0,1]". Then, 



(a) T{xi,. . . ,x„) = oo <=^ YJ^=i l{i^>>o} < & 

1 

i=l 

2. Let (yi,...,2/„) 6 (0,oo]". Then, 

(a) Tiyi, y„) > ^{yi=oo} < b; 



(h) Setting y, = Tixk ■ k^i), we have V ^{vi<oo} = blr 



(b) Setting x[ = Tiyt -k^i), we have ^ l{x'<oo} = U{x\,. ..,x'^)-b^^ l{j,^=oo}- 

Proof. Properties l.a and 2. a are straightforward from the definition of F and F. Regarding property 
1.6, set K ^ ^ {i ^ [n] : > 0} . If the sum is non-zero then there must be an i satisfying both yi > 
and Xi < oo. By l.a, this implies K > b. Conversely, if 7^ > 6 then < oo for every i £ [n]. We have 
just shown 

T!. n 



_ ^r^J2\I\=b.l3iYlkeI^k 

- ^{K>b} 2^^^ yf — 

= &l{i^>6}, 

where the second equality is obtained by replacing = F(a;fc : fc 7^ i) by its definition. For property 
2.6, set L = ^ {i G [n] : yi — 00} . When L = 0, 2.6 boils down to formula (fTTj). The case 1 < i < 6 can 
then be obtained by successively setting each of the L variables to 00, which amounts to condition on 
the presence of the corresponding ground elements. For L > b, both sides of the equation are zero. □ 

Lemma 7 Assume that P ^ Q ^ P' . Set s ^ Pii'^Y), t = Q{{oo}) and s' = P'({0}=). Then, 

1. s UtU s' ; 

2. P' <P ^ MiP) < His) ; 

3. P' >P=^ MiP) > His). 

In particular, if P ^ (6 o e)(P) then MiP) ^ His) and W is) = 0. 
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Proof. In the whole proof, N denotes a generic random integer with law tt, iV a generic random inte- 
ger with law TT, X, Xi, X2, ■ ■ ■ generic (0, 00]— valued random variables with law P, Y, Yi, 1^2, • ■ • generic 
[0, 1]— valued random variables with law Q, and X' , . . . generic (0, 00]— valued random variables 

with law P' . We use the convention that all variables appearing under the same expectation are inde- 
pendent. With these notations, properties l.a and 2. a in Lemma [6] give 



N 



N 



and 



'(X'>0)-P^l{y, = «,}<6 



i=l 



which, in view of the definition of /, yields exactly the first claim of the Lemma. Now, using property 
l.b and 2.6, we respectively obtain the two following identities : 



E 



XY 



1 + XY 



^{Y<oo} 



riGN 



XT {Xi, . . . , Xn^l) ^ ^ 
l+XT{Xi,...,Xn^l) {r(Xi,...,X„_i)<oo} 



E/ N „ XT {Xi, . . . ,Xn~l) 
Tr[n)nE 



neN 



\' ~ ±T • • ' -J IL~ ± / -. 

1 + XT{X,, X„_i) Mr(^i.-.^"-i)<-} 



X,T{Xk:k^i) ^ _ 



N 



El 



{X,>0} 



> h 



h{\-g{s)). 



E 



X'Y 
1 + X'Y 



L{y<oo} 



7r(n)nE 

new 

E '^w^ 



yr(yi,...,y„,i) . 
i + yr (yi,...,y„_i)- 

i + rr (yi,...,F„_i 



St 



- 1 + y.r (Ffc : A; 7^ i) 



-{F<oo} 

;l{y<oo} 
1-{yi<oo} 



E[c/(ri,...,rjv)] -E 



AT 



00} 



2M{P)-b{l-g{s))-csfis). 



Since the mapping is increasing in x, claims 2 and 3 follow. □ 

Proof of TheoremlB Fix s e [0, 1] satisfying H'{s) = 0, i.e. (/ o /)(s) = s. Define G P([0, 1]) to be 
the Bernoulli distribution with parameter s, and set then iteratively 

F,:+i-(eoe)(p„^) 

for all neN. By part 1 in Lemma [71 Pf is a distribution on [0, 1] satisfying PidO}"^) = s. Since Pq 
is the largest such distribution, we have Pf < Pq. But both Q and O are decreasing, so by immediate 
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induction, the sequence {P^JneN is non-increasing. Thus, the hmit = hni„^oo | Pn exists in ^([0, 1]). 
Setting Soo = P^{{0}'^), we claim that 

1. P'* is a fixed point of o 8 ; 

2. A/(P-^) = i7(soo); 

Part 1 follows from the continuity of 8 and 8. Part 2 is guaranteed by Lemma [T) Part 3 is a consequence 
of the fact that < Pg . Finally, for each n e N, we have Pni{OY) = s and P^ > P^+i, so Lemma [3 
guarantees that M{P^) < H{s). Letting n — > oo yields exactly part 4. 

We are now in position to prove the equivalence between the historical minima of H and the solutions 
to (8 o 8)(P) — P. If s is a historical minimum, then parts 3 and 4 force Soo = s so P = P'* is a 
solution satisfying M{P) = H{s), as desired. Conversely, we must show that any solution P is in fact 
of the form P'* for some historical minimum s. Set s = P({0}^), which satisfies H'{s) = by Lemma 
[T] We have P < P^ for any n G N : this holds for n = because PJ is the largest element of P([0, 1]) 
such that PodO}"^) = s, and it then inductively extends to all n S N by monotony of 8 o 8. Letting 
n ^ oo, we obtain P < P" ; but by Lemma[7]we also have M(P) = H{s) = M{P''). Thus, P ^ P" (Af 
is decreasing). Finally, if i < s is any historical minimum then clearly Pq < PJ, which implies P^ < P'^ . 
In fact the inequality is strict, because P*{{QY) = t < s = P''{{0}'^). Applying the decreasing mapping 
M yields H(t) > H{s), which shows that s is a historical minimum. 
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